Role of initial quantum correlation in transient linear response 
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The linear transient response of a two-level system coupled with an environmental system is stud- 
ied under correlated and factorized initial conditions. We find that the transient response in these 
cases differs significantly from each other, especially for strong system-environment interaction at 
intermediate temperatures. This means that it is necessary to pay attention to the initial conditions 
chosen when analyzing experiments on transient linear response, because the conventional factor- 
ized initial condition results in an incorrect response, in which the quantum correlation between the 
relevant system and the environmental system is disregarded. 
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The linear response to an external field has been in- 
tensively studied to provide information on the charac- 
teristics of condensed matter. Since the early work on 
nuclear magnetic resonance (NMR) by Wangsness and 
Bloch we often use a factorized initial condition of 
ps ® Pr, where the statistical operator ps for the rel- 
evant spin system, and that pr for the environmental 
system separately stay in each equilibrium state. This 
means that we conventionally neglect the quantum cor- 
relation between the system and the environment at the 
initial time which is justified only for a sufficiently weak 
system-environment interaction. 

We should pay attention to the fact that such factor- 
ization is sensitively reflected mainly in the transient be- 
havior of the linear response to an external field, because 
the system feels the above mentioned difference in the 
initial condition after a sudden switch-on of the external 
field. However, its effects have been discussed from var- 
ious viewpoints such as (1) the reduced dynamics of an 
open quantum system I2mj0j, (2) the complete positivity 
of dynamical maps |lll4l5| , and (3) the nonlinear optical 
response [l^]. The effect of a correlated initial condition 
on the linear response has been investigated by focusing 
only on the stationary susceptibility [l7|- Because the de- 
velopment of experimental techniques enables us to study 
transient behavior in a short time region, it might be im- 
portant to investigate how the factorization of the initial 
condition affects the transient response under a suddenly 
applied external field. In this paper, we answer this ques- 
tion using an exactly solvable model, i.e., a two- level sys- 
tem that linearly interacts with an environmental system 
(bosonic thermal bath). 

The Hamiltonian of the matter system consists of a 
relevant two-level system, a bosonic bath, and the linear 
interaction between them, which are given in the form 



with 



Us = S 1 |1)(1| + S |0)(0|, 
Kr = y^TiLo k b\b k , 

k 

n SR = 5><? fc (4 + &fc)|i)<i|, 



(2) 



where -Eo(i) is t ne energy of the lower (upper) state of 
the relevant system, is the frequency of the bosonic 
bath mode, bl and are its creation and annihilation 
operators, and is the coupling strength between the 
excited state and the boson mode. We study the tran- 
sient behavior of a two-level system abruptly excited by 
an external field. The matter-field interaction Hamilto- 
nian is given by 

Hp(t) = -i/x-Ee(t)|l)(0|e-^ t -i/x*-S^)|0)(l|e i ^ i , 

(3) 

where p is the transition dipole moment, E is the am- 
plitude of the external field, uj p is the frequency of the 
external field, and 8(t) is the step function. It is con- 
venient to use a canonical transformation in terms of 
S = exp[S|l)(l|] with B = J2k(9k/u k )(b k - b{), because 
this allows us to eliminate the system-bath interaction 
from the matter Hamiltonian in the form 



E[\l){l\ 



(4) 
(5) 



with E[ = Ei — fij^kidk/^k)- The matter-field interac- 
tion is transformed as 



W P {t) = S^n P (t)S 



= --p-E9(t)\l){0\e 



(6) 



H = Hs + Hr + Hsr, 



(1) 



Under the application of an external field, the time evo- 
lution of the density operator of the total system, p(t) 
for t > 0, is given by 

p{t) = Sd n,t U(t)p'(0)UHt)e^ n,t S\ (7) 
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where we define p'(0) = S^p(0)S and 

U(t)=T + exp[-- H P '(t')dt'] 



(8) 



with H P '(t) = e^'^pWe*^''. Now we consider the 
linear response by making an approximation as 



U(t) w 1 - - / Hp (t')dt' 



(9) 

The initial statistical operator in the equilibrium state 
is transformed as 

p'(0)=p' 1 |l)(l|+^|0)(0|, (10) 

which gives the transformed density operator in the form 

p'{t) = &p(t)S 
«e^'Vi|l)<l|+^|0><0| 



, 1 



M-^(G(t'H-PiG(O)|l)(0| 



- 1 



-p*.E(GHt')p[-p' GHt'))\0)(l\}}e-^ t . 

(11) 

where we define G(t) = e st(t) with B*(t) = 
e ^'*Bt e -H-«'t and Aw = {E' 1 -E )/h-uj p . The induced 
dipole moment under the application of an external field, 



/J(i)=Tr[(/x|l)(0|+/r|0}(l|)p(t)] 



(12) 



is evaluated by using Eq. (fTTjl for each initial condition 
with trace operation for total system Tr which consists 
of trace operation over system (environmental) variables 
Tr S (fl)- 

As an initial condition of the total system, we consider 
the following two cases: (1) the total system is in an 
equilibrium state and (2) the two-level system and the 
environmental system are in separate equilibrium states. 

Case (1): Correlated initial condition 
The initial state p(0) = exp[— f¥H\ = ptot is trans- 
formed as p'(0) = exp[— P'H' S ] exp[— PHr] with Z = 
Tr[exp[— /3H]], which means that 



-PE'i 



Pi 



PR 



Po 



e- pE °p R 
7' ' 



(13) 



where we defined Z' s = Trg[exp[— j3"H' s ]) and p R = 
^exp[-/3H fl ] with Z R = Tr R exp[-f3H R }. 

Case (2): Factorized initial condition 
For the factorized initial state p(0) = ps ® Pr, with 
Ps = -^exp[-/3Hs] and Z s = Trs exp[-fiHs], the 
transformation gives 



Pi 



e~^G{0)p R G\Q) 
Zs 



Zs 



Po 



e- 0Eo p R 



Zs 



(14) 



with p' R = J?- exp[— H' R ] and Z' R = TrR_exp[— H' R \ , where 
we define 



(15) 



This means that for the factorized initial condition, the 
excited state of the two-level system effectively interacts 
with a displaced bosonic reservoir. It is noted that the 
initial statistical operator for case (1) (Eq. ([13])) has a 
factorized form, and that for case (2) (Eq. (fH)) ) has a 
correlated form. This is simply because of the canonical 
transformation, and should not be confused with corre- 
lated and factorized conditions in the original represen- 
tation. It is also noted that we have another choice of 
the initial decoupled state as p(0) — Trg[p tot ] Cg>TrR,[ptot]> 
which will be discussed subsequently. 

Defining p,r m \(t) as an induced dipole moment for case 

(m) with m = 1,2 apart from a factor of 2p(fl* ■ E)/K, 
we obtain , 

M(m)(f) = \A( m )(t)\ COs(uJ p t - <j>f m )(t)), (m = 1, 2) 

(16) 

where A( m ^(t) is defined by the following equations, 
i rt 

= / *V Aw(t - t) {e- /3B i*i(t-t') 
z s Jo 

-e- ?Eo ^l{t - 1')}, (17) 

A m {t) = -L I dt'e iAu{ t-^{e-P B ^ 2 {t,t') 
^s Jo 

-e- 0Eo ^l(t-t')}, (18) 

and 4>( m )(t) is its argument. With a coupling spectral 
function as h(oj) = ^Zufft^i^ — lu^), the renormalized 
energy is given by E[ = E\ — U J °° dujh{uj)/w. ^\{t) and 
*2 (*,*') in Eq. (JT7]) and Eq. JTSJ are obtained as 



#i(t) = exp[-£(t)-i 



h(u) 



sin(art)], (19) 



with£(t) = r o °° duj(h(u})/uj 2 )(l + 2n{uj))(l~cos{ujt)) and 
* 2 (M') = Vi(t-t') x 



exp[ 



duj 



h(u) 



2i(shi(u)t ) — sin(cjt))]. 



(20) 



As is seen from Eq. ()13|) ~ ()20p , the difference between 
the two cases in the transformed representation results 
in (1) a shift in the transition frequency ujq, and (2) the 
displacement of the boson operator, which appears in the 
difference between and ^ 2 {t,t') in Eq. (JT5JI and 

In the following, we show the numerical evaluations of 
the transient response for the Ohmic coupling spectral 
density, h(u>) = su>e~ u / Uc . In this case, the renormalized 
energy is given by E[ — E\ — Tislu c , and we find that the 



3 



imaginary part of the exponent of Eq. (|T9| and (1201) is 
integrated with a formula as j Q du](h(uj)/uj 2 ) sin(cji) = 
s arctan (w c . This means that the difference between 
the time evolution of ^i(t) and ^(i, i') is written as 
exp[— 2is(arctan (u> c t') — arctan (uj c t))]. 

In Fig.[TJ we show the time evolution of the intensity of 
the dipole moment under the application of an external 
held for fe^T = IQhujQ, s = 1, u> c = ujq/5, and uj p = ujq. 
Time is scaled as t = w^t. We hnd that the evolution 
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FIG. 1: (color online) Time evolution of the induced dipole 
moment for ksT — 10/iajo, s = 1, D c = Wo/5, and uj p = Do- 
Dark gray (blue) and light gray (red) lines represent cases 
(1) and (2), respectively. Time is scaled as t — Uot. (a) 
time evolution of intensity of dipole moment |/i( m )(f)| 2 , (b) 
the time evolution of amplitude |.A( m )(i)| and (c) the time 
evolution of phase, (j>( m )(t) with m = 1,2. 

of case (1) approaches stationary oscillation much faster 
than that of case (2). The characteristic time of the 
slower rise for case (2) depends on u> c , which comes from 
the difference between *&x(t — t') and ^2(t,t'), i.e., the 
arctan part. We show in Fig. [ljb) the time evolution of 
amplitude, |.A( m )(t)| with m — 1,2, where we find the 
amplitude for case (2) approaches larger value than that 
for case (1). The overestimation of the dipole moment for 
the factorized initial condition results from the fact that 
the decoherence due to quantum correlation is not taken 
into account in case (2). As shown in Fig. [He), we find 
that the two cases have different phases, <p/ m \(t) with 
m = 1,2, for their oscillations. The feature originates 
from the effective displacement for the factorized initial 
condition shown in Eq. (fl4l) . 

When the temperature is decreased to ksT = hujQ, we 
find that the difference in the amplitudes of the station- 
ary oscillations for the two cases increase (Fig. [Ha) and 
(b)). And the asymptotic values of phase 0( m )(t) differ 
each other between m = 1,2. 

However, as shown in Fig. [3] for a lower temperature 
hujo/5, we find that differences in amplitudes and phases 
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FIG. 2: (color online) Time evolution of the induced dipole 
moment for a lower temperature case ksT = hujQ. Other 
parameters and evaluated quantities are the same as in Fig. [TJ 



of the two cases become very small. This is because at 
sufficiently low temperatures, the initial occupation prob- 
ability for the upper level is very small, and the initial 
quantum correlation is disregarded. 
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FIG. 3: (color online) Time evolution of the induced dipole 
moment for ksT = Hljo/5. Other parameters and evaluated 
quantities are the same as in Fig. [T] 

From these considerations, the effect arising from the 
initial correlation becomes significant in the intermediate 
temperature region, where the energy difference between 
the two levels becomes comparable to thermal excitation. 

In order to see the difference in the transient response 
under the two types of initial conditions at the high tern- 



4 



perature limit, we set n(uj) w 1/((3Huj) 
analytic form of £(t) in Eq. P^|) as 



which gives an 



arctan(u; c i)- 
for the Ohmic spectral function. 



(2-/3^ c )log(l+w c 2 i 2 ) 



(21) 



In Fig. [4] (a), we show the time evolution of intensity of 
dipole moment \^{ m ){t)\ 2 for fc^T = 10 4 7iu;o with ui c — 
wo/50, s = 1, and w p = wq- Fig. H] (b) and (c) show the 
time evolution of amplitude |A( m )(i)| and phase 4>( m )(t) 
with m = 1,2, respectively for the same parameters as 
in Fig. Ufa). 
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FIG. 4: (color online) Time evolution of the induced dipole 
moment for a high temperature case, k B T = 10 4 hu} , with 
lo c = Wo/50, s = l, and lo p = wo- The evaluated quantities 
are the same as in Fig. [T] 



We clearly find that a slow rise for case (2) for a smaller 
lo c . Generally, l/u) c corresponds to the correlation time 
of the system-environment interaction, and the slow rise 
time determined by the long correlation time is an arti- 
fact arising from neglecting the initial quantum correla- 
tion. In contrast to the former cases of lower tempera- 
tures, we find that the asymptotic values of A/ m \(t) and 
<A(m)(£) f° r m — 1 , 2 agree with each other very well in 
the long-time region. 



In summary, we compared the transient linear response 
under two types of initial conditions, i.e., with and with- 
out initial correlation. We found that the response for 
these cases differed from each other for strong system- 
environment interaction at intermediate temperatures. 
Although the present analysis was limited to a two-level 
system linearly coupled with a bosonic environmental 
system, the initial quantum correlation may generally al- 
ter the linear responses of various types of materials. 
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